Kelvin-wave cascade and dissipation in low-temperature superfluids vortices 
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We study the statistical properties of the Kelvin waves propagating along quantized superfluid 
vortices driven by the Gross-Pitaevskii equation. No artificial forcing or dissipation is added. Vortex 
positions are accurately tracked. This procedure directly allows us to obtain the Kevin-waves 
occupation-number spectrum. Numerical data obtained from long time integration and ensemble- 
average over initial conditions supports the spectrum proposed in [L'vov and Nazarenko, JETP Lett 
91, 428 (2010)]. Kelvin wave modes in the inertial range are found to be Gaussian as expected by 
weak-turbulence predictions. Finally the dissipative range of the Kelvin-wave spectrum is studied. 
Strong non-Gaussian fluctuations are observed in this range. 
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Superfluid turbulence has been the subject of many 
experimental and theoretical works for the last decades. 
In particular a lot of progress has been done in experi- 
mental techniques and it is now possible to realize tur- 
bulent Bose-Einstein condensates (BEC) jl , turbulent 
flow with 3 He [5J |3] and visualize vortex filaments in 4 He 
[3] . As in classical hydrodynamic turbulence [5] , a turbu- 
lent Kolmogorov cascade has been observed experimen- 
tally and numerically. In superfluids, this takes place 
at scales larger than the mean inter-vortex distance I 
[SHE]. At very low temperature, when damping due to 
mutual friction is negligible, it is believed that dissipa- 
tion at small scales is carried by phonon radiation which 
dissipates energy into heat [9]. At scales smaller than i 
the energy is transferred down by a series of reconnec- 
tion processes of quantized vortices that excite waves on 
the filaments. These perturbations called Kelvin waves 
(KW), have been known for more than one century [10] 
in the broader context of fluid dynamics. These waves 
obey a set of non-linear equations where the energy trans- 
fer towards small scales is carried by a wave-turbulence 
cascade. How the energy is distributed along different 
scales is crucial for the understanding of the dissipative 
processes in superfluids. The energy spectrum of such a 
cascade is not yet fully-determined, except in the limit 
of small amplitude KW, where the theory of weak turbu- 
lence is applicable [TT] . However, a heated debate on the 
locality of KW energy transfer has taken place in the last 
years [HHTTj- Two different groups, Kozik & Stvistunov 
[18] and L'vov & Nazarenko [19], starting from the very 
same equations and by using the same theory, have de- 
rived two different spectra (hereafter KS and LN spectra 
repectivelly). The origin of this controversy is mainly due 
to a symmetry argument by KS (tilt of a vortex line) that 
eventually leads to a vanishing vertex in the perturbative 
expansion. This leads to locality in the energy transfer 
and makes the 6-wave interaction theory realizable. The 
energy spectrum found by KS is 

£ KS (fc)~ e 1/5 K 7/5 fc- 7/5 , (1) 



where e is the energy flux, n the circulation quantum and 
k the wavevector. This symmetry argument was ques- 
tioned by LN and they claimed that the energy transfer 
is non-local. They derived an effective 4-wavc interaction 
theory that leads to the energy spectrum 

E LN (k) ~ ^e 1 / 3 *- 2 / 3 ^ 5 / 3 , (2) 

where ~ (1/ K ) J Ejj$(Je)dk is the mean-square angular 
deviation of the vortex from its straight-line configura- 
tion. More technical details on the controversy can be 
found in pMCF]. The exponent 7/5 = 1.4 and 5/3 w 1.67 
of and ^ , are respectively supposed to be universal, 
but their relatively close values makes it difficult to nu- 
merically elucidate which theory is correct. A number 
of numerical works supporting both theories have been 
published in the last years but none presenting strong 
arguments to settle this controversy [13 [2Ql [21] . These 
works are all done in the framework of the vortex filament 
(or equations derived from them) with an ad-hoc dissipa- 
tive mechanism. It is worth mentioning that in the case 
of strong wave-turbulence, when the local slope of KW 
is order one, weak-turbulence breaks down and a differ- 
ent spectrum was proposed by Vinen ct al. 22 . The 
Vinen spectrum Ey[ nen (k) ~ & > was derived under the 
assumption of critical balance where the linear and non- 
linear time-scales of KW are of the same order. Finally, 
It was suggested by E. Sonin [16] that no universally can 
be expected for the KW spectrum. 

In this Letter, we address the small-amplitude KW 
cascade problem by performing direct numerical simula- 
tions of the Gross-Pitaevskii equation (GPE). The GPE 
formally describes the dynamics of a weakly-interacting 
BEC at very low temperature. It is also expected to at 
least qualitatively reproduce the dynamics of superfluid 
Helium vortices. As the Gross-Pitaevskii (GP) vortices 
can naturally radiate and excite phonons no artificial dis- 
sipation is needed, making it the natural framework for 
studying the KW cascade and the low-temperature dissi- 
pative processes. We first show that vortices coexist with 
small-scale thermalized phonon-waves. Then, the (II?) 
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KW occupation-number spectrum is precisely obtained 
and data is found to support the wave-turbulence pre- 
diction (LN) [19]. The KW spectrum is analyzed within 
the dissipative range and an exponential decay is found. 
Finally, the probability distribution function (PDF) of 
KW amplitudes is observed to be Gaussian in the iner- 
tial range in contrast with the power-law tails observed 
for modes in the dissipative range. 

The GPE describing a homogeneous BEC of volume V 
with (complex) wave-function ip is given by 

*4^-^vV + .9HV, (3) 

where m is the mass of the condensed particles 
and g = Airah 2 /m, with a the s-wave scatter- 
ing length. Madchmg's transformation ip(x,t) = 

exp [z^0(x, t)] relates the wave-function ip to a 
supcrfluid of density p(x, t) and velocity v = V</>, where 
k = h/m is the Onsager-Feynman quantum of velocity 
circulation around the tp — vortex lines. Equation ^ 
conserves the total energy H = J(|^|VV>| 2 + §|^| 4 )(2x 
and the total number of particles N — J |-0| 2 c?x. When 
Eq. (pi is linearized around a constant ip = ipo, the sound 

~ 1/2 

velocity is given by c = (g\ipo\ 2 /m) with dispersive ef- 
fects taking place at length scales smaller than the coher- 

1/2 

ence length £ = (h 2 /2m\ipo\ 2 g) that also corresponds 
to the vortex core size. In the GPE numerical simulations 
presented in this Letter the density p = mN/V is fixed 
to 1 and the physical constants in Eq. ^ are determined 
by the values of £ and c — 2. The quantum of circula- 
tion h/m has the value c^/\f2. Numerical integration of 
Eq.([3| is performed by using a standard pseudo-spectral 
code with an exponential time-splitting temporal scheme 
in a cubic box of length V 1 ^ 3 = 2n. Resolutions used are 
listed in Table |T] below, the value of £ is chosen as small 
as possible but well resolving vortex lines, and ensemble- 
averaging is done over 30 initial conditions to reduce fluc- 
tuations. 

To address the KW problem, we use the simplest con- 
figuration allowed by periodicity consisting of an array 
of four alternate-sign vortices. To obtain a clean initial 
condition and reduce the initial phonon emission, in a 
first step, an exact stationary solution of the GPE with 
straight vortices is numerically obtained by a Newton 
method [23] . Vortices are separated by a distance tt and 
can be considered isolated when £ — > 0, as the resolution 
is increased. Then, a KW is introduced slightly perturb- 
ing the vortex. The initial KW reads 

x{z) = A cos (nz + <j>l) , y(z) = A ^ cos (nz + <p v n ), 

n n 

(4) 

where <j)\? ,y ^ are random phases and n = 2, 3. A 3D vi- 
sualization of the the density p(x) is displayed on FigJT^ 
at t = 10. The KW is clearly observed with the (red) iso- 
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FIG. 1: a) (color online) 3D visualization of the density \%jj\ 2 in 
the sub-box [0, ir} 2 x [0, 2ir}. In red an isosurface of the vortex 
and a (orange) density plot shows sound waves, b) Temporal 
evolution of total, incompressible kinetic, compressible kinetic 
and wave energy, b.l) and b.2) zoom of incompressible and 
compressible kinetic energy respectively, c) Incompressible ki- 
netic, compressible kinetic and wave energy spectra. Dashed 
lines display k 2 and k~ 5 ^ 3 power-law scalings. 



surface. Phonon-waves correspond to the (orange) cloud 
that is a density-plot of p(x) in a narrow threshold cen- 
tered around the mean density value p = I. To quantify 
the vortical and wave energy of the configuration, we use 
the standard hydrodynamic energy decomposition, ob- 
tained by using the Madelung transformation (see [23] 
for details). The total energy is thus decomposed in two 
terms: the incompressible kinetic energy i?£ in containing 
the contribution of vortical structures and the energy of 
phonon-waves E wliv = ££ in + E [nt + E q , where E q , E iat , 
-^kin' are t ne quantum, internal and compressible kinetic 
energy respectively. Figure [T]b displays the temporal 
evolution of ££ in , E^ n , E wav and E tot = El in +E wav . Ob- 
serve in Fig{l]3.1 and Fig(T|3.2 that their temporal evolu- 
tion rapidly reach a (quasi-)statistical stationary regime. 
The same energy decomposition can be applied to the 
energy spectra that are displayed on Fig[lJ: at t = 10. 
The energy spectrum of the compressible kinetic energy 
presents at large wavenumbers a fc 2 -equipartition regime. 
This range is also present in the initial condition albeit 
with smaller values. It rapidly reaches the stationary 
state shown in Figjlfc showing that thermalized waves 
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coexist with vortices. We note that E^ in <C E wav , and 
hence the large-scale GPE dynamics is mainly driven by 
vortices setting a clean configuration to study the KW 
cascade problem. 

The energy spectra displayed on Fig{T]c present a scal- 
ing compatible with fc -5 / 3 , however it cannot to be asso- 
ciated with Komogorov turbulence as the scale separation 
V 1 l' i > £ is not realized (here V 1 ^ ~ £ for all times). 
The scaling could be explained by the presence of a KW 
cascade and predictions ([l} or pj), as the principal con- 
tribution to energy of the fluid (see FigjlJ is mainly due 
to vortices. However the relationship between the purely 
ID (D is the dimension) KW spectrum and 3D (hydrody- 
namical) energy spectra is not clearly stablished. To ex- 
plicitly study the KW cascade, we numerically track the 
coordinates (x(z),y(z)) of the vortex. For each value of z 
the equation ijj(x(z),y(z)) = is solved by using a New- 
ton method with the wave- function tfj obtained from ^ . 
Derivatives of the fields at inter-mesh points are obtained 
by Fourier interpolation. This allows us to accurately ob- 
tain the vortex coordinates with a precision much larger 
than the one given by the mesh size or any other mesh- 
interpolation. Once the coordinates are obtained, it is 
possible to compute the KW vortex length, mean curva- 
ture and in particular the (ID) KW occupation- number 
spectrum (hereafter KW spectrum) defined by 



n(k) = \w{k)\ 2 + \w(-k)\ 2 , 



(5) 



where w(k) is the Fourier transform of w(z) = x(z) + 
iy(z). The KW spectrum allows us to directly construct 
the KW energy Ekw and the energy dissipation e: 



E K w = y^^(fc)n fc 



dE 



KW 



dt 



(6) 



where ui(k) is the KW dispersion relation. The KW 
dispersion relation can be approximated by uj(k) = 
C '(k I 'Air)k 2 , where C is numerical constant, which even- 
tually depends logarithmically on As the previ- 
ous quantities are computed from a ID signal, they 
all present strong fluctuations. To reduce fluctua- 
tions, ensemble average is performed over the phases 
of the initial condition Q. Figures [2^, and [2)3 show 
the temporal evolution of the total vortex length L = 
j v/l-i- \d z w(z)\ 2 dz and the mean curvature K — 
J \d z w(z) x d zz w(z)\/\d z w(z)\ 3 dz, normalized by their 
initial values. Note that both quantities slightly fluctu- 
ate with time. Finally, on Fig(2j) the temporal evolution 
of the KW energy is displayed. Note that the energy 
fluctuates, especially after the arrival of phonon waves 
coming from neighboring cells near t ~ tt/c ps 1.5. The 
solid (red) line presents the energy averaged over tem- 
poral windows of width At = 2; the decrease in energy 
is apparent. In the inset of Figj2]D, the temporal evolu- 
tion of the energy dissipation e is shown. Note that the 
energy dissipation fluctuates and shows some negatives 
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FIG. 2: (color online) a) Vortex length L(t) (blue circles) and 
curvature K(t) (red crosses) normalized by L(0) = 6.51 and 
K(0) = 0.44. b) KW energy E KW {t) (circles). The (red) 
solid line displays i?Kw(£) averaged over temporal windows 
of width At = 2. Inset: energy dissipation e. 



values as well. This can be related to the presences of 
phonon waves that excite KW at small scales. However 
the temporal average is positive as more energy is radi- 
ated than absorbed by the vortex. The value of e is not 
precisely determined here, as it depends on the constant 
C and on the way temporal averaging is done, but in any 
case its mean value remains positive. 

We now turn to the KW spectrum. Two kinds of sim- 
ulations are presented: The first, trying to enhance the 
scale separation between I/ 1 / 3 and £ and thus obtaining 
a larger inertial range. The second, concerns the dissi- 
pative range of the KW spectrum and then presenting a 
large number of modes between £ and smallest resolved 
scale V^ 3 /N z . The different runs and their details are 
listed on Table [Q 

Let us focus now on the inertial range of the KW cas- 
cade. The two KS and LN wave-turbulence predictions 
for the KW spectrum read 



«ks(k) = ^Yy/g , n L N(K) 



47rC LN e 1/3 

^2/3^.11/3 ' 



(7) 



where Cks and Cln are numerical constant. 

The temporal evolution of the KW spectrum is dis- 
played on Fig|3]a for the Run III. KW are rapidly exited, 
and populate all wavenumbers. Energy arrives to scales 
small enough to be efficiently dissipated creating a steep 
decay zone usually called dissipative range in hydrody- 
namic turbulence [5j . As dissipation by phonon emission 
is very weak [25], the system reaches a (quasi)-stationary 
state where a clear inertial range can be observed. To 
reduce fluctuations in the spectrum, time average dur- 
ing the stationary regime has also been performed. Fully 
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FIG. 3: a) Temporal evolution of KW spectrum, Run III. b) 
Time averaged KW spectra, RunS I-V. Dashed line displays 
the power-law fits. 



averaged spectra of Runs I-IV are displayed on Fig|3]3. 
An inertial range with a power-law scaling is clearly ap- 
preciated for almost one decade. The exponent to, ob- 
tained from a fit k~ rn , is shown on Table []] with their 
respective errors. Note that Run II and III have been 



Run 


N ± 


N z 


TV 

1 * rca 


n 


4 


A 


m 


I 


256 


128 


31 


3 


0.025 




3.85 ±0.24 


II 


256 


128 


31 


2 


0.025 


4£ 


3.682 ±0.13 


III 


256 


256 


11 


2 


0.025 


4£ 


3.753 ±0.17 


IV 


512 


256 


1 


2 


0.0125 


4£ 


4.116 ±0.56 


V 


128 


512 


11 


2 


0.1 


4£ 





TABLE I: List of runs. N± and N z are the resolutions in the 
perpendicular and parallel directions respect to the vortex. 
TVrea is the number of realizations, n corresponds to the num- 
ber of initial KW modes and m corresponds to the exponent 
k~ m of the KW spectrum. 

performed with the same physical parameters but dif- 
ferent resolutions in the z direction resulting in almost 
identical spectra. The choice of resolution used in Run 
II allows us to speed up computations, and thus to get 
better statistics, at the expense of a reduced dissipative 
range. For all runs the exponent is slightly larger than 
the one predicted by the two weak-turbulence results Q 
and presents a variation of a 5%. However, for all runs 
data supports the exponent —5/3 — 2 predicted by LN, 
that it is within errors, and excludes the —7/5 — 2 KS 
prediction. Experimentally, one usually have access only 
to the (3D) kinetic energy, that for small amplitude KW 
is dominated by the singularity of the velocity at the 
vortex core. A singularity cannot transfer energy and 
the KW cascade is thus crucial for understanding low- 
temperature dissipative-mechanisms of supcrfluids. The 
progress that is being achieved in visualization techniques 
of vortex filaments [3] , will maybe allow in the future to 
experimentally addressed this issue. 



We now turn to the dissipative range of the KW spec- 
trum, that takes place at scales much smaller than l/£. 
For such small scales it is known that dispersive ef- 
fects of phonon waves slowdown the dynamics produc- 
ing a bottleneck and quasi-thermalization |26j of modes 
in an intermediate range. This effect was numerically 
observed for large values of £fc ma x, where fc max is the 
largest wavenumber. A natural question is if this slow- 
down can affect the dissipative process of the KW cas- 
cade? If excitations of high wavenumber phonons are 
difficult, one could expect that dissipation of KW by 
sound emission should be reduced at such scales. Fur- 
thermore, at such small scales KW amplitudes should 
be somehow decoupled from the large-scale movement of 
the vortex and must be in direct interaction with out of 
equilibrium phonon waves. To investigate such a configu- 
ration we have performed simulations with a large value 
of £&max = 17 (Run V). For such a configuration, the 
inertial-range of the KW cascade is not clearly identi- 
fied in the KW spectrum presented in Figj4^L. A similar 




FIG. 4: a) Temporal evolution of KW spectrum for the 
Run V. Inset: temporal evolution of the exponential decay 
rate S(t). b) PDF of KW amplitudes in the inertial range 
20 < k < 40 (Run III), c) PDF of KW amplitudes in the 
dissipative range 80 < k < 100 (Run III) and 30 < k < 170 
(Run V). Inset: same PDF in Log-Log, the power-law to~ 7//3 
is drawn for reference 



behavior to the one observed in the fully 31? dispersive 
bottleneck is observed at very early times. KW stop to be 
populated at wavenumbers larger than fcj = 27r/£ ~ 60. 
An equipartition regime of KW numbers for k < and 
faster than exponential decay for k > k^ is observed. The 
similarity with the 3D dispersive bottleneck must be con- 
sidered carefully. First, the dispersion relation of KW 
is dispersive at all wavenumber unlike the phonon dis- 
persion relation (or Bogoliubov) where dispersive effects 
become important at scales smaller than £. Secondly, in 
the 3D dispersive bottleneck an equipartition of energy 
is observed and not equipartition of phonon occupation 
number. Furthermore, the statistics of the KW ampli- 
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tudes in this range (30 < k < 60, t — 0.1) is not Gaus- 
sian as expected in thermal equilibrium. At later times, 
the equipartition range is destroyed and a large-fc expo- 
nential decay is stablished showing the analyticity of the 
fields. The rate S(t) (obtained with a fit n k ~ e -2 *^*) 
is displayed in the inset of Fig|4]a. It fluctuates in time 
with a minimum value <5 m i n = 0.007. 

Finally we study the statistics of KW amplitudes in the 
inertial and dissipative range. The PDF of KW ampli- 
tudes varying at scales inside the inertial range for Run 
III is obtained by filtering in Fourier space and keeping 
modes between the range 20 < k < 40. The normalized 
PDF is displayed in Fig|4jD together with a Gaussian dis- 
tribution. The quasi-Gaussian behaviors (except for the 
very long tails) is manifest as expected from weak tur- 
bulence predictions. The square amplitude |u>(z)| 2 con- 
sequently presents exponential tails. In the dissipative 
range (for 80 < k < 100) the PDF shows a strong non- 
Gaussian character as apparent on Fig|4jc. Better statis- 
tics are obtained when considering Run V and modes 
between 30 < k < 170. The PDF has power-law tails 
as shown on the inset of Figj4jj. Also, as is usual for 
turbulent flows 5J, an asymmetry is found when com- 
puting the skewness (u> 3 ) / (w 2 ) 2 / 3 = —0.15. Recently, in 
Biot-Savart simulations [27 , a crossover between Gaus- 
sian and non-Gaussian statistics was found in the veloc- 
ity field taking place at the scale I. Here, for KW, the 
crossover takes place at the scale £ <C t. At scales smaller 
than £ the KW are somehow decoupled of the large scale 
dynamics prescribed by wave-turbulence and are in direct 
interaction with a strongly fluctuating superfluid velocity 
field, as the one found in the GP simulations of Ref. [35J . 

The behavior of KW at very small scales is an impor- 
tant issue, especially at very low temperature where dissi- 
pation by mutual friction is absent. In all vortex filament 
models, some small-scale artificial dissipation is needed 
to avoid energy pile-up at the smallest resolved scale. 
Although vortex-filament models are not concerned with 
such smalls scales, how the energy is dissipated in those 
models can affect the vortex dynamics and it would be 
important to check if the dissipative mechanisms used 
are consistent with the dissipation produced by phonon 
radiation in the GP results presented on Figj4] 

To conclude, the KW cascade has been studied by 
performing Gross-Pitaevskii simulations of a slightly- 
perturbed straight vortex. Precise tracking of vortex 
coordinates and ensemble averaging allow us to the de- 
termination of a clean KW spectrum that supports the 
wave-turbulence prediction of L'vov and Nazarenko (LN) 
[T5] , At scales smaller than the vortex size, an exponen- 
tial decay of the KW spectrum was found. The PDF of 
KW was found to be quasi-Gaussian in the inertial range 
and to have power-law tails in the dissipative range. Pre- 
cise tracking of GP vortices can be used in the future to 
improve the understanding of vortex thermal- waves inter- 
action by using, for instance, the truncated (projected) 



GPE [53], where mutual friction effects are present. 
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